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Abstract 

The possibility of a traumatic bone fracture in space is a concern due to the observed decrease in 
astronaut bone mineral density (BMD) during spaceflight and because of the physical demands of the mission. 
The Bone Fracture Risk Module (BFxRM) was developed to quantify the probability of fracture at the femoral 
neck and lumbar spine during space exploration missions. The BFxRM is scenario-based, providing predictions 
for specific activities or events during a particular space mission. The key elements of the BFxRM are the 
mission parameters, the biomechanical loading models, the bone loss and fracture models and the incidence 
rate of the activity or event. Uncertainties in the model parameters arise due to variations within the population 
and unknowns associated with the effects of the space environment. Consequently, parameter distributions 
were used in Monte Carlo simulations to obtain an estimate of fracture probability under real mission 
scenarios. The model predicts an increase in the probability of fracture as the mission length increases and 
fracture is more likely in the higher gravitational field of Mars than on the moon. The resulting probability 
predictions and sensitivity analyses of the BFxRM can be used as an engineering tool for mission operation and 
resource planning in order to mitigate the risk of bone fracture in space. 
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Topics to cover 



Overview of Integrated Medical Model (IMM) and 
Bone Fracture Risk Module (BFxRM) 

Definition of Fracture Risk Index (FRI) 

Library of biomechanical models used to estimate 
load on bones during activities and events 

Decrease of Bone Mineral Density (BMD) in space 
and relationship between BMD and ultimate strength 
of bone 

Model uncertainty 

Earth based validations of models 

Sample results - probability of fracture on moon and 
Mars missions 

Conclusions 
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Integrated Medical Model (IMM) 
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The Integrated Medical Model (IMM) is a tool for quantifying the 
probability and consequences of medical risks 

Integrate best evidence in a quantifiable assessment of risk 

Identify medical resources such as skills, equipment, and supplies 
necessary to optimize mitigation strategies. 
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Bone Fracture Risk Module (BFxRM) 
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Fracture Risk Index (FRI) 

• The ratio of the load experienced by the bone during an activity to the 
ultimate strength of the bone. 

- An FRI of less than one indicates that the bone should be strong enough 
to support the load 




- An FRI of greater than one indicates that there is a significant risk of bone 
fracture. 

Loads experienced by the bone are The ultimate strength of bone is found from 
estimated with biomechanical models testing the strength of cadaver bone 



Chaffin DB, Baker WH, A biomechanical model for analysis of symmetric 
sagittal plane lifting, AIIE Transactions , 2(1), March 1970, pp. 16-27 




Yoganandan N, Pintar F, Sances A, Maiman D, Mykelbust J, Biomechanical investigations of the 
human thoracolumbar spine, In Biomechanics of Impact Injuries and Injury Tolerances of the 
Abdomen, Lumbar Spine and Pelvis Complex , Edited by SH Backaitis, Society of Automotive 
Engineers, Inc., 1995, pp. 97-114. 
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Library of biomechanical loading models 



Femoral Neck - Fall to the side 


Hip mass 

Stiffness and damping 
of hip pad and ground 



S. N. Robinovitch, W. C. Hayes, and T. A. McMahon, ’’Prediction of femoral impact 
forces in falls on the hip,’’ J. Biomech. Eng , vol. 113, no. 4, pp. 366-374, Nov. 1991. 


Lumbar Spine - Trunk flexed, 
holding a load 



Lumbar Spine - Fall, 
landing on two feet 


Upper body 
mass 

Stiffness and damping 
of lumbar spine 

Pelvis and 
leg mass 

Stiffness 
of leg 


Foot mass 

Stiffness and damping 
of ground 



A. Schultz, G. B. Andersson, R. Ortengren, R. Bjork, and M. Nordin, ’’Analysis and 
quantitative myoelectric measurements of loads on the lumbar spine when holding 
weights in standing postures,” Spine, vol. 7, no. 4, pp. 390-397, July 1982. 


K. J. Chi and D. Schmitt, ’’Mechanical energy and effective foot mass 
during impact loading of walking and running," J. Biomech., vol. 38, 
no. 7, pp. 1387-1395, July2005. 
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BMD Loss in space over time 




BMD DoE = BMD value on the day of the event 

BMD Start = BMD at the beginning of the mission 

BMD Loss = The amount of BMD loss prior to the day of the event 


Data used to determine slope includes LSHA Data and 
Published Data and takes into consideration uncertainty 
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Relationship between BMD and Ultimate Load 
of bone for different loading conditions 




BMD (g/cm 2 ) 


K. Singer, S. Edmondston, R. Day, P. Breidahl, and R. Price, "Prediction of thoracic and lumbar vertebral body compressive strength - Correlations with Bone Mineral 
Density and vertebral region," Bone, vol. 17, no. 2, pp. 167-174, 1995. 


E. N. Ebbesen, J. S. Thomsen, H. Beck-Nielsen, H. J. Nepper-Rasmussen, and L. Mosekilde, "Lumbar vertebral body compressive strength evaluated by dual-energy X- 
ray absorptiometry, quantitative computed tomography, and ashing," Bone, vol. 25, no. 6, pp. 713-724, Dec. 1999. 


D. P. Lindsey, M. J. Kim, M. Hannibal, and T. F. Alamin, "The monotonic and fatigue properties of osteoporotic thoracic vertebral bodies, "Spine, vol. 30, no. 6, pp. 
645-649, Mar. 2 00 5. 


B. S. Myers, K. B. Arbogast, B. Lobaugh, K. D. Harper, W. J. Richardson, and M. K. Drezner, "Improved assessment of lumbar vertebral body strength using supine 
lateral dual-energy x-ray absorptiometry," J. Bone Miner. Res., vol. 9, no. 5, pp. 687-693, May 1994. 
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BFxRM uncertainty 

• Monte Carlo and Latin Hypercube simulations performed to determine 
most likely probability since: 

- The system parameters (i.e. astronaut mass, reference BMD level, BMD 
loss per day, ultimate strength/BMD, anthropometric values, physiological 
stiffness and damping constants, angle of trunk flexion, load lifted, etc.) are 
defined as distributions over a range of values. 

- The event could happen on any day during the mission 




Normal Distribution 




Damping constants 




Damping constant of lumbar spine (Ns/m) 
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Earth based validations- Static lumbar spine 

model 



Comparison of lumbar spine loading 
calculations 



Comparison of FRI calculations 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 

Condition 

Condition 1 are young subjects 
Condition 2 are elderly subjects 


□ Duan et al. (2001) 

□ Calculated 


Y. Duan, E. Seeman, and C. H. Turner, "The biomechanical basis of vertebral body fragility in men and women," J. Bone Miner. Res., 
vol. 16, no. 12, pp. 2276-2283, Dec. 2001. 


M. L. Bouxsein, L. J. Melton, III, B. L. Riggs, J. Muller, E. J. Atkinson, A. L. Oberg, R. A. Robb, J. J. Camp, P. A. Rouleau, C. H. 
McCollough, and S. Khosla, "Age- and sex-specific differences in the factor of risk for vertebral fracture: a population-based study using 
QCT," J. Bone Miner. Res., vol. 21, no. 9, pp. 1475-1482, Sept.2006. 
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Earth based validations- Static lumbar spine 

model 



Comparison of Ultimate Load 
vs. Age 


Comparison of % FRI above 
1 vs. Age 




M. L. Bouxsein, L. J. Melton, III, B. L. Riggs, J. Muller, E. J. Atkinson, A. L. Oberg, R. A. Robb, J. J. Camp, P. A. Rouleau, C. 
H. McCollough, and S. Khosla, "Age- and sex-specific differences in the factor of risk for vertebral fracture: a population- 
based study using QCT, " J. Bone Miner. Res., vol. 21, no. 9, pp. 1475-1482, Sept.2006. 


M. Biggeman, D. Hilweg, S. Seidel, M. Horst, and P. Brinckmann, ‘‘Risk of vertebral insufficiency fractures in relation to 
compressive strength predicted by quantitative computed tomography, ’’Euro J Rad, vol. 13, pp. 6-10, 1991. 
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Earth based validations - Dynamic lumbar 

spine model 

Comparison of Ground Reaction Force calculations 




* 

* 

* 


Seegmiller et al. (2003) 

Arampatizis et al. (2001 ,2001 ,2004) 
McNair et al. (1999) 

Kwok et al. (2003) 

Calculated mean 
Mean + 3 Std 
Mean - 3 Std 


Height (cm) 

J. G. Seegmiller and S. T. McCaw, ’’Ground Reaction Forces Among Gymnasts and Recreational Athletes in Drop Landings,” J. Athl. Train., vol. 38, no. 4, pp. 311- 
314, Dec. 2003. 

A. Arampatzis, G. P. Bruggemann, and G. M. Klapsing, ’’Leg stiffness and mechanical energetic processes during jumping on a sprung surface," Med. Sci. Sports 
Exerc., vol. 33, no. 6, pp. 923-931, June2001. 

A. Arampatzis, F. Schade, M. Walsh, and G. P. Bruggemann, "Influence of leg stiffness and its effect on myodynamic jumping performance," J. Electromyogr. 
Kinesiol., vol. 11, no. 5, pp. 355-364, Oct.2001. 

A. Arampatzis, S. Stafilidis, G. Morey-Klapsing, and G. P. Bruggemann, "Interaction of the human body and surfaces of different stiffness during drop jumps," Med. 
Sci. Sports Exerc., vol. 36, no. 3, pp. 451-459, Mar. 2004. 

P. J. McNair and H. Prapavessis, "Normative data of vertical ground reaction forces during landing from a jump," J. Sci. Med. Sport, vol. 2, no. 1, pp. 86-88, 

Mar. 1999. 

P. Kwok, W. Kong, K. Kasturi, C. Lee, J. Hamill, “A biomechanical study on the parachute landing fall, ” 17th AIAA Aerodynamic Decelerator Systems Technology 
Conference and Seminar, 19-22 May 2003, Monterey, CA., AIAA 2003-2149. 
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m b k Earth based validations - Dynamic lumbar 
~ spine model 


Comparison of fracture prediction for a fall height distribution 




Our simulations predicted an FRI above 1 for 34.2% of the trials. Goonetilleke 
found 29.7% of falls in his study resulted in fracture. 


U. K. Goonetilleke, "Injuries caused by falls from heights," Med. Sci. Law, vol. 20, no. 4, pp. 262-275, Oct. 1980. 
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Example results 

Probability of fracture of the lumbar spine by a male or female astronaut 
due to lifting a load with the trunk flexed during an EVA during a long 
duration, Lunar mission. 

Male Female 





Gender 

Mean Probability 

Standard 

Deviation 

5% 

95% 

Male 

3.19e-4 

1.17e-4 

1.84e-4 

5.36e-4 

Female 

3.28e-4 

1 .36e-4 

1 .8e-4 

5.85e-4 
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Example results 

Probability of fracture of the lumbar spine by a male or female astronaut 
due to a 1m fall during an EVA during a long duration, Martian mission. 



Male Female 


1 00,000 T rials Frequency View 97,927 Displayed 



1 00,000 T rials Frequency View 97,360 Displayed 
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Gender 

Mean 

Probability 

Standard 

Deviation 

5% 

95% 

Male 

2.64e-3 

5.36e-3 

5.54e-5 

1 . 1 9e-2 

Female 

3.02e-3 

6.00e-3 

5.97e-5 

1.39e-2 
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Conclusions 

• A model has been developed that bounds the 
uncertainty associated with the risk of bone fracture in 
space. 

- Integrative approach accounting for extenuating factors 

• Equipment and Vehicle 

• Bone Health 

• Training and Operations 

• The model can be used to predict the most likely 
probability of bone fracture in space. 

- “what if scenarios 

• What if reduced gravity is osteo-protective? 

• What if the FFD is reduced to t-score of -1 .25? 




• The model can be used as a useful engineering tool 
during mission planning. 
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Future Work 


Wrist fracture risk assessment 

Renal stone formation risk assessment 

Insomnia and circadian rhythm upset risk 
assessment 
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Calculating Bone Ultimate Structural Strength 

Ultimate Structural Load 
Capacity for Loading Conditions 


Posterolaterial fall: 
UL Reduced up to 
0.8% per Degree 




Apply UL attenuation for 
load direction 


Use BMD correlations to 
Estimate UL 


Based on appropriate 
ex vivo test data 


Estimate Time Course to and Degree 
Of Bone Loss at Skeletal Location 
On day of loading 


State of Bone at 1g 
Pre-Flight DEXA-BMD 




Linear or Exponential Model 


Maximum Loss Est. 
With Pop. Variability 


NHANES DATA - Represents Pre- 
Flight Bone Health, FFD Standards 
And Reference Max BMD Condition 
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Calculating Loading in Reduced Gravity 
Suit Environment 



Mass & Padd ing 


Resultant Skeletal Load 





ssi* Determine Load Additive or 

Active Response Attenuation Factors 



Scale Load to Gravity 
Using Appropriate Methods 


Estimate of Load 
Biomechanics 
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Uses the change in momentum 
Includes additional mass 


Loading Event Occurs 
From Specified 
Activity or Incident 



UPDATE: Newest models use 
simulations based on Robinovitch 
type loading simulations 




Represents a perceived 
loading state during on 
surface activities 


www.nasa.gov 21 





Probability of 
FN Fracture 
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Apollo 
Data 

Fall Rate: 0.35/hr and a = 0.066 
Pr(Postlat): 0.0517 and a = 0.0404 


Probability bone 
will fail to 
support load 


FRI Estimates 
From BFRM 
Bone Loss 
Bone Strength / 
Quality 

Loading Levels in 

Hypo-g 

Mission 

Characteristics 

Equipment / Suit 

Characteristics 


Published 
Data 
Relating 
FRI and 
Fracture 
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Estimated upper and lower bounds: 
FRI To Probability of Fracture 
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® Schaffner 
© Results 


“Smell” Test Validation 



£0.000 Trials Frequency View 19.956 Displayed 



Pre-flight estimate of FRI for Unhindered Posteriolateral Fall 
i.e. a fall to the side and slightly backward 
Male in 1 g with ~1m fall heights 
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Probability of Fracture Due to Side Falls 
Male on Extra Vehicular Activity 



Data Shown for Mars: 540D Surface Mission 



25 : 000 Trials 


Frequency View 


23 : 669 Displayed 


M:L 



4.00000E-3 6.00000E-3 

Probability of Fracture 
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T1 
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-Q 
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£ 
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0 


► 1 1 .63200E-6 


Certainty: |90.000 


4 |l 09094E-2 


Mission 

Fracture Probability 

Std 

5th Percentile 

95th Percentile 

Lunar: 8D Surface 

1.50E-4 

1.15E-3 

3.30E-07 

5.36E-04 

Lunar: 170D Surface 

1.94E-4 

1.54E-3 

3.47E-07 

6.15E-04 

Mars: 40D Surface 

1.44E-3 

7.66E-3 

1.15E-06 

4.85E-03 

Mars: 540D Surface 

2.47E-3 

9.95E-3 

1.68E-06 

1.15E-02 


Lateral/Posteriolateral Fall heights range from .25m to ~1m 
Bone loss not attenuated by partial gravity 
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Model Sensitivity 



Lunar: Long 

1 00,000 T rials Contribution to Variance View 


Sensitivity: Trochanter - Unhindered fall 

-11.0% 0.0% 11.0% 22.0% 33.0% 44.0% 



Mars: Long 

1 00 : 000 T rials Contribution to Variance View 


Sensitivity: Trochanter - Unhindered fall 

-10.0% 0.0% 10.0% 20.0% 30.0% 



• The suit attenuation characteristics and the impulse scaling factors 
produce the most sensitivity 

• Interesting to note that 

- Successful reaction to the fall is the next most driving factor 

- Bone loss rates are not as significant for lunar missions 

- Reference BMD produces more sensitivity to the calculation than rate of 
bone loss in both scenarios 
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